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LONG-TERM  GOALS 

The  long-range  goals  of  this  project  are  to  develop  an  improved  turbulence  closure  model  that  takes 
explicit  account  of  surface  wave  effects  and  to  quantify  the  effect  of  Stokes  drift  on  the  surface  current 
signature  of  High  Frequency  (HF)  radar  systems. 

OBJECTIVES 

The  objectives  of  this  project  are  to:  1)  deploy  a  surface  buoy  with  turbulent  flux  and  surface  wave 
measurement  capabilities  within  the  footprint  of  the  HF  radar  network  in  Monterey  Bay  and  2) 
formulate  a  new  turbulence  closure  scheme  based  on  surface  wave  dynamics  that  can  be  evaluated  using 
data  from  the  flux  buoy. 

APPROACH 

This  project  is  extending  the  measurement  and  modeling  efforts  already  underway  in  the  Monterey  Bay 
region  to  include  the  direct  measurement  and  modeling  of  surface  waves.  It  is  combining  three  efforts  to 
assemble  a  better  surface  flux  parameterization:  1)  the  Ly  et  al.  wave  modeling,  2)  the  Davidson  et  al. 
meteorological  flux  observations,  and  3)  the  NOPP/ICON  coastal  modeling. 

WORK  COMPLETED 

Three  separate  deployment  of  the  NPS  Flux  Buoy  have  taken  place  for  30-60  day  periods:  Fall  1999, 
Winter  2000,  and  Fall  2000  (co-sponsored  by  an  ICON  follow-on).  The  formulation  of  a  new 
turbulence  closure  scheme  based  on  surface  wave  dynamics  has  been  completed  and  the  initial 
implementation  of  the  scheme  has  been  attempted  within  the  POM-based  models  of  Ly  (NAval 
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postgraduate  school  ocean  Model;  NAM)  and  Shulman  (Innovative  Coastal-Ocean  Observing  Network 
model;  ICON). 

RESULTS 

The  instrumented  NFS  FLUX  buoy  was  assembled,  checked  out  and  deployed  in  Monterey  Bay  to 
provide  in  situ  measurements  of  surface  waves  and  atmospheric  forcing  by  K.  Davidson  of  the  NFS 
Meteorology  Department  (see:  http://www.met.nps.navy.mil/~jonesk/BuoyHistory/ppframe.htmI.  In 
addition  to  bulk  wind  and  temperature  observations,  critical  turbulent  flux  data  on  wind  stress  and 
directional  surface  wave  spectra  were  collected  at  the  mooring  site  within  Monterey  Bay.  Of  particular 
interest  here  is  the  directional  wave  spectra,  which  are  being  used  in  the  model  to  test  the  wave  age 
parameterizations.  These  continuous,  2-D  wave  data  are  also  being  used  to  assess  the  impact  of  Stokes 
currents  on  high  frequency  (HF)  radar  backscatter  measurements.  A  sample  wave  spectrum  from  the 
Fall  1999  deployment  is  shown  in  Figure  I. 
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Figure  1.  Directional  wave  spectrum  for  1400  UTC  on  5  October  1999. 

A  new  turbulence  closure  with  surface  wave  parameterization  (Ly,  1995;  Ly  and  Garwood,  2000;  Ly  et 
al.  2000a)  was  developed  for  the  oceanic  and  atmospheric  models  taking  surface  ocean  waves  into 
account.  This  turbulence  closure  has  been  developed  based  on  a  new  concept  of  turbulent  energy 
sources  related  to  both  shear  production  (classical  theory)  and  breaking  waves  using  the  semi-empirical 
turbulence  and  similarity  theories,  and  a  linear  theory  of  breaking  waves.  The  turbulence  scheme  with  a 
new  wave  parameterization  has  been  tested  using  the  NFS  I-D  air- wave-sea  interaction  coupled  model 
and  new  data  of  dissipation  distribution  in  the  ocean,  and  other  air-sea  interaction  characteristics,  such 
as  wave-dependent  drag  coefficient,  roughness  length,  and  momentum  fluxes.  The  numerical  studies  and 
data-model  comparisons  show  that  the  new  theory  can  predict  (better  than  the  classical  theory)  wave- 
dependent  turbulent  kinetic  energy  dissipation  (Figure  2)  in  the  ocean  (and  it  is  expected  in  the  marine 


atmospheric  turbulent  layer).  The  wave-dependent  drag  coefficient  (Figure  3)  and  wave-dependent 
roughness  length  also  agrees  well  with  the  observed  data  (Ly  and  Garwood,  2000).  The  most  important 
point  is  that  a  close  comparison  of  the  new  theory  with  all  available  data  was  made  and  the  first  full 
application  of  the  model  results  to  ocean  circulation  models,  such  as  the  Princeton  Ocean  Model 
(POM)  and  NAM  (Ly  et  ah,  2000a),  was  conducted,  beginning  the  move  toward  a  true  coupled 
circulation-wave  system. 


Figure  2.  Various  measurements  (symbols)  of  nondimensional  dissipation,  epsilon,  as  a  function  of 
nondimensional  depth.  Data  points  adapted  from  Agrawal  et  al.  (1992)  along  with  the  classic  wall- 
layer  theory  distribution  and  a  range  of  results  from  the  new  wave  parameterization. 


Figure  3.  Dependence  of  model  drag  coefficients  on  10-m  height  winds  and  wave  ages  equal  to  10, 
20,  and  30.  Three  other  curves  with  "very  young”,  "mature”,  and  "full"  from  observational  data  by 

Smith  etal.  (1992). 

The  new  turbulence  closure  with  a  full  wave  parameterization  has  been  used  to  develop  the  NAM 
ocean  circulation  model  with  multi-block  grid  features,  which  allow  better  handling  of  complicated 
bottom  topography,  coastlines  and  multi-scale  physical  phenomena  in  modeling  domains.  Multi-block 
grid  techniques  also  have  superior  properties  in  parallel  computing  in  comparison  with  the  traditional 
single  block  grid  and  domain-decomposition  techniques  (Luong  et  ah,  2000abc).  The  model  has  also 
been  used  to  develope  a  barotropic  tidal  model  for  the  Monterey  Bay  region.  These  results  have  been 
written  up  as  an  invited  paper  for  the  International  Hydraulic  Engineering  Software  VIII  Conference, 
2000  (Ly,  et  ah,  2000b). 

As  part  of  this  project,  a  modification  of  the  new  surface  wave  parameterization  for  idealized  wave 
height  and  wave  age  fields  has  been  also  incorporated  within  the  ICON  modeling  and  data  assimilation 
effort.  Our  first  version  of  a  new  PROFQ  subroutine,  along  with  some  modifications  to  the  MAIN 
program,  have  been  transferred  to  Igor  Shulman  and  James  Lewis  of  that  project.  Also,  a  new  version 
of  subroutine  PROFU  with  an  equation  of  energy  fluxes  for  taking  drift  current-wave  relationship  is 
being  transferred  to  the  ICON  modeling  group  for  testing.  It  is  noted  that  the  wave  parameterization 
was  partly  modified  for  use  in  POM  and  other  models  utilizing  mixing  length  theory.  In  this  way. 


POM  can  take  the  wave-turbulenee  relationship  (without  turbulenee  dissipation)  into  aeeount.  (At  this 
time,  only  the  NPS  NAM  version  ean  take  a  full  wave-turbulenee  relationship  into  aeeount.) 

The  NAM  box  model  was  used  in  a  eireulation-wave  eoupling  study  with  idealized  wave  height  and 
wave  age  fields  (Ly,  et  ah,  2000a).  This  study  foeused  on  the  sensitivity  of  the  eurrent  field  to  the 
surfaee  waves.  It  also  demonstrates  the  eapability  of  NAM  in  reprodueing  an  upwelling  feature  for  the 
idealized  California  eoastal  region.  Simulation  showed  that  surfaee  eurrents  are  strongly  sensitive  to 
surfaee  waves  in  a  eomplex  wind-wave-turbulenee-eurrent  relationship.  Simulations  also  show  that 
oeean  models  that  do  not  take  surfaee  waves  into  aeeount  ean  overestimate  eurrents  up  to  50  %. 

IMPACT/APPLICATIONS 

The  likely  impaets  of  this  researeh  will  inelude  improved  understanding  of  the  effeets  of  surfaee  waves 
on  the  measurement  of  surfaee  eurrents  by  HF  radar  installations  along  with  an  expanded  eapability  to 
inelude  the  effeets  of  surfaee  waves  within  turbulenee  elosure  sehemes  used  by  primitive  equation 
models.  The  applieation  of  these  results  should  lead  to  fundamentally  improved  models  and  more 
aeeurate  remotely  sensed  surfaee  eurrents. 

TRANSITIONS 

The  transition  opportunities  are  related  to  next-generation  oeean  eireulation  models  and  eoupled 
atmosphere-wave-oeean  models.  These  models  will  require  improved  parameterizations  of  the  effeets 
of  surfaee  waves  on  turbulent  transport  near  the  oeean  surfaee.  In  addition,  the  information  about  the 
effeets  of  surfaee  waves  on  HF  radar  measurements  will  be  needed  to  eorreet  observed  surfaee  eurrent 
fields  prior  to  assimilation  into  oeean  eireulation  models. 

RELATED  PROJECTS 

This  projeet  represents  and  extension  of  the  efforts  underway  as  part  of  the  Monterey  Bay  eomponent 
of  the  National  Oeean  Partnership  Program  (NOPP/ICON).  It  also  eomplements  the  ICON 
eoordination  with  the  field  tests  of  the  ONR-sponsored  Autonomous  Oeean  Sampling  Network 
(AOSN)  eondueted  in  Monterey  Bay  in  August  2000.  Information  about  this  projeet  and  the  broader 
ICON  and  AOSN  efforts  ean  be  found  at:  http://www.oc.nps.navy.mil/~icon  . 
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